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Abstract 


When an energetic particle (kinetic energy >05 MeV) originating from a radio- 
active deca\ or a cosnic ra' traverses the active regions of semiconductor devices 
used in integrated circuit (IC) chips, it leaves along its track a high density 
electron-hole plasma The subsequent decay of this plasma bv drift and diffusion 
leads to charge collection at the electrodes large enough in most cases to engender 
a false reading, hence the name single-event upset (SEU) The problem of SElVs is 
particularly severe within the narsh environment of Jupiter's radiation belcs and 
constitutes therefore a matter of concern for the Calileo mission 

In the following \e shall anal\/e the physics of an SEU event in some detail Cur 
mam conclusion is this Cuing to the predominance of nonlinear c nace charge 
effects and the fact that positive (holes) and negative (electrons) charges must 
be treated on an equil footing, analytical models for the ion 1 zed-charge collection 
and their corresponding currents as a function of time prove to be inadequate 
even in the simplest case of uniformly doped, abrupt p-n junctions in a one-di- 
mcnsional geometry The necessity for full-fledged computer simulation of the 
pertinent equations governing the electron-hole plasma therefore becomes -»mpera- 
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1 Introduction It is well known that an energetic paiticle (electron, 

a-particle, cosmic ray, etc ) produces electron-hole pairs when traversing a 
semiconductor The number of pairs produced per cm of track length ** s given by 

"o • - 4k a> 

where -dE/dx is the stopping power or linear energy transfer (LET) of the particle 
and w the mean energy spent to produce one electron hole pair For Si at 300°K 
this number is w = 3 6 eV independent of the type of particle^ For example, a 

Q 

5 MeV ct-particle produces Nq - A v 10 electron-hole pairs per cm, a number which 

2 ) 

can easily be calculated using standard stopping power tables The track of 

the energetic particle or, more specific, the volume in which elec t ron-hole pairs 

are generated can well be approximated by a cylinder of length R, where R is the 

*) 

range of the fast particle possessing a circular base area corresponding to a * 

-5 3) 

radius of 0 1 vm = 10 cm . In the above mentioned example, this gives for the 
number density of electron-hole pairs generated by the fast a-particle 

N = N /it x 10~ L0 = 1.25 x 10 18 cm" 3 (2) 

a number considerably larger than the usual doping levels of extrinsic material 

16 — 3 ** 

(N a - 10 cm ). Since dE/dx of eq* (1) depends on x, will in general also 

be a function of x, a fact which has been analyzed in more detail by Yaney et al 

for a-particles^ . 

Here is not the place to delve into this particular matter any further The fact 
that energetic particle beams generate copious amounts of electron-hole pairs m 
semiconductors is being used extensively m high energy particle detectors^ and 


If the range is larger than the semiconductor sample size, the latter must of 


course be used 
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m materials research (determination of diffusion lengths by means of the scanning 
electron microscope for instance^) The property of energetic particles to 
readily ionize semiconductor material used so advantageously m the examples just 
given becomes a detrirent in mary other cases For instance RAM cells of modern 
computers consist of a number of transistors and diodes which m turn are fabri- 
cated from Si p-n junctions Information is encoded as different voltages or 
different charge states at appropriate nodes of the circuitry These voltages 
and charge states are constantly changing as information is processed It has 
been realized for some time that an energetic particle crossing a RAM cell may 
upset the circuitry by virtue of the generation of electron-hole pairs which when 
collected at sensitive nodes may give rise to false readings^ Since then a 

fairly large body of literature has arisen dealing with this problem which is now 
generally known under the name of "single-event upset" (SEU)^ 

On the next pages we shall formulate the problem, show it^ inherent difficulties 
and show that analytical models are inadequate Th^s state of affairs is caused 
by an enhanced collection of carriers at the nodes via drift in the originally 
space charge neutral region of a p-n junction This is referred to in the liter- 
ature as "funnelling H Incidentally, funnelling, which is a typical high level 
injection effect, was discovered by computer modelling rather than analytical 
methods and cane as a surprise to the community 


2. Formulation of the problem We consider for simplicity here a p-n junction 

+ 20 -3 
consisting of a narrow N region (donor concentration ^ 10" cm ) joined to a 

15 - 3 

wider P region (acceptor concentration = IQ cm ) Again for simplicity we 

\ 

take as the energetic particle a 5 MeV a-particlv, ' wmch penetrates the planar 



^5 MEV a-particles were chosen because both analytical and computer models 
exist for this case (see below) , 
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device normally A 5 MeV n-parricle possesses a range of 24 um in Si 

-12 

is 0 = v/c = 0 052 A device 20 urn thick will be traversed in 1 28x10 sec 

During this time electron-hole pairs will be created along the track in the 
— 8 — 1 5 

amount of Nq » 4 \ JO cm At most Q = 8 x 10 elementary charges can be 

collected from this ion track at any node For static RAM's, both the total charge 
and the rise time (the current dQ/dt) are important Consequently, the time 
history of the excess carriers must be known For times much greater than 
thermalizat ion times '«10 ^ sec) carrier transport proceeds via drift and 
diffusion The election and hole currents are then given by 


3 = q u nE + q D Vn (3a) 

in ^ n ~ 1 n 

and 

Jp = q Up pE - q VP (3b) 


-19 

Here q is the elementary charge 1 6 y 10 A see, v and D are the mobility and 

n n 

diffusion constant for electrons, a and D those for holes, E is the electric 

P P 

field The number densities of electrons n and of holes p satisfy the continuity 

8 ) 

equations 
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Here G and R signify total generation and recombination rates In Si, recombination 

proceeds predominantly via recombination centers with energy levels within the 

9) 

forbidden gap According to the Shockley-Read-Hall theory , G-R, both for 

electrons and holes G - R and G - R respectively, become functions of the 

n n p p 1 J 

density of charged recombination centers m such a manner that^^ 
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R + G - G 
n n p 


(5) 


3 


yt~ r — * 


' ' - "j 
OF PuC»< Q j / 


G f " ' f ' 


In cq (5) t. and G are the generation rates of electrons and holes associated 
n p 

with recombination centers whereas the difference 


G - C = C - G 
n n p p 


( 6 ) 


signifies generation rates by otter mechanisms (for instance fast a-particles) 

where electrons and holes are generated at an equal rate n* and n^ are the 

number densities of positively and negatively charged recombination centers 

Io cLose the s>stems of equations (3) to (6), the electric field L must be 

related self-consistently to n and p For time-varying currents and electron- 

—12 

hole concentrat ions with characteristic time constants ^ 10 1 sec, the electric 
field is appt oximately irrotational and we can put 


E - - 


(7) 


where is an electrostatic potential, a function of both position and time The 
final connection is made through Poisson's equation 


^ (n ' P + N a - N d + n T - n T > 


( 8 ) 


which constitutes a relationship between the number density of charged carriers, 

charged acceptors and donors, charged recombination centers, and the gradient 

of the electric field, and is nothing more than a manifestation of Gauss's law"^ . 
-14 

r q = 8 854 v io A sec/V cm is the dielectric constant of vacuum, e = 10 8 for 
Si signifies the relative dielectric constant of the medium. and are 

assumed to be known time-independent functiois of position. At room temperature 
and for non-degenerate material they are equal to the total accepter and donor 
Np densities respectively From eqs (4), (5), (7) and (8) it follows that the 
total current density j is conserved 
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where 


J = 0 


j = j + j -f,c L 

i in ip o ^ 


(9a) 


(9b) 


The set of equations (3) to (9) forms the most general basis f rom which to 
attack the SEU problem However, before this can be even contemplated , boundary 
conditions and initial conditions must be formulated \s we have seen earlier, 
the i-particle generates electron-hole pans in a narrow coLumn rather suddenlx, 
l e , in a time short compared to current rise times due to drift ind diffusion 
of the excess carriers The appropriate initial condition nay then be formulated 
m the following way Initiallv, at t me t = 0, the electric field distribution 
and the electron as well as the hole concentrations are those prevailing before 
the ^-particle strike except for an additional density of electron-hole pairs 
within a right circular cvhnder along the t-particlc tiack with a radius of about 
0 1 jin and a magnitude for the excess carrier density computed from eq (I) 
Although the formulation of an appropriate initial condition is rjLhcr straight- 
forward as we nave seen, boundar'y conditions are less easily accessible Ibis 
is because, as we shall see m the next section, we are dealing here with 
phenomena which have an appreciable content of amplitudes in the 10^ to 10^ Hz 
frequency range of their Fourier spectrum and it is not at all clear that 
boundary conditions applicable to static or low frequencx situations are valid 

in the high frequency regime The static boundary conditions for the ohmic con- 

12 ) 

tacts of a transistor have however been recently applied to a c signal 

13) 

analyses with seeming success The static boundary conditions at ohmic 

i u 12 > 

contacts assume equilibrium 


np = n 


( 10 ) 


where n is the intrinsic carrier concentration, and charge neutrality 
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Furthermore the line integral of the electric field between any two contacts must 

*) 

be equal to the potential difference e>isting between the contacts 

2 

E dS = V ^ - V (12) 

1 

At free surfaces of cue semiconductor the usual boundary conditions in the steads 
state are 


in *12 = ‘ S(n • n 0 } 


ip 212 " " S(p “ p 0 } 


(13) 


with n^^ the outward normal and S the surface recombination velocity and 

are the equilibrium values for electron and hole concentrations respectively 

14 ) 

Now, according to Many et al., S is a function of the number density of recom- 
bination centers at the surface and m a time-dependent situation, part l culir 1 \ 
for high frequencies, these states will empty and fill according to certain rate 
equations governed among other things by the time dependence of electron and hole 
densities at the surface* - ^ As a result, the common surface recombination 
velocity S of eqs. (13) splits into two different time-dependent expressions 

S (t) and S (t) valid for electrons and holes separately Iortunatol\ v Si sur- 
n P 

faces used m computer hardware ire passivated and their surface recombination 
velocity is small. It is possible to use as a first approximation 


in 


il2 


H\2 ~ ® 
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Fhe line integral (12) must include the built-in potentials of any -junction 


crossed by the path of integration 
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or in other words, the normal component of electron and hole currents vanishes 
at free surfaces Eqs (10), (11), (12) and (14) constitute the necessary 
boundarv conditions to implement a solution of the lasic equations (3) to (9), 
given the initial condition discussed in the text 

3 Methods of solution . The set of equations and boundary conditions, eqs (3) 

to (14), constitutes indeed a formidable non-linear system which cannot be 
solved in full generality, even with sophisticated computers Therefore approxi- 
mations are inevitable In order to find sensible approximations we must look 
it the physical system described by the equations m more detail Confining 
ourselves just to one p-n junction in reverse bias as an element of a static 
RAM cell, we know that within the depletion layer approximation, the junction 
consisting of a narrow n + laver (the emitter) and a wide, moderately doped p-type 
bulk (the base), admits of three distinct regions The quasi-neutral emittet, 
a depletion layer and Emails a neutral (field free) base region r or a one- 
dimensLonal abrupt junction, the depletion layer width is gi\en by 


W 



(V, 4- \ ) 


i /: 


(15) 


whore V is the built-in voltage and V the reverse bias voltage is m turn 

given bv 


V 


b 


kT 

q 


n(N A yn;) 


(16) 


The depletion la\er contains a strong elettric field v» mb is due to the polarized 
dipole la\er If an x-coord mate is introduced positive in the direction from 
the n + -t\pc to the p-t>pe material, the electric field may be expressed as 


E = E 

m 


(1 


(17) 
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M , 1/2 


(V + V) J 

D 


Eqs (17) and (18) ire again valid in the depletion la>er approximation. For 

15 -*3 

t>[ical Si junctions, assuming an acceptor density of N =10 cm , a donor 
19 -3 

density of = 10 cm , room temperature T = 300 C K and a reverse bias V = 

8 volt, we find for W and E 


W = 3 24 am, E = 5 43 x 10 V/cm 
m 


As we see, the electric field drops linearly from a rather high value to zero 


at the edge of the depletion layer 


\ 5 MeV a-particle striking the p-n junction at normal incidence will traverse the 
device and Lea^e behind a column of a high density electron-hole Dlasma as already 
discussed Pait of this plasma finds itself in the depletion layer where it is sub- 
jected to a large electric field, part of it has been generated within the field- 
f**ee bulk of the device The generation of free carriers happens almost instantan- 
eously, so that initially the field configuration remains what it was before the 
jL-stnke. But immediately following the u-strike that part of the plasma generated 
by the u-particle which occupies the depletion layer will be torn apart by the strong 
electric field and drift toward the electrodes, while the other part within the base 
will diffuse toward the electrodes The charges will be collected ay the electrodes \ 
and it can be seen that there are two mechanisms, which are represented by two tune 
constants, governing the flow of carriers drift and diffusion. The mobility of 
holes gp is about twice as small as the mobility of electrons u n in Si. Therefore 
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M will determine Che time constant for prompt collection Since the speed of 

P 

holes is equal to j E, we have as an estimate for the time constant governing cLift 

t s I'/u E = l(f 11 sec (20) 

t pm 

using the values for W and given befoic (eq (19)) ard using for the hole 

2 15 -3 

mobility u = 600 cm /V sec valid for the doping density of h = 10 cm Die 
P A 

diffusion exhibits another time constant which is considerable larger and ma> 
be estimated b> 


S 4 D i\\2 
P 


2 x 10 sec 


where we used for the width of the base 1 = 10 M m and for the diffusion constant 

2 

for holes D = 15 5 cm /sec \ ] e see therefore, that charge collection proceeds 
P 

fast by drift (prompt charges) and slowly by diffusion (delayed charges) This 

qualitative picture has been borne out both by experiment and computer simulat ion"^ 

Because of the anticipated small time constant (20) for prompt charge collection, 

the computer program^^ could be simplified drastically by neglecting the recom- 

-k-k ) 

bmat lon-generat ion terms G - R of eqs (4; In this case numerical solutions 

of the pertinent equations of the previous section become possible, but require 
extensive use of computer time During the course of these calculations the effect 
of "field funnelling" was discovered -\s the electron-hole pairs situated along the 
a-track m the depletion layer are separated by the strong electric field (18), 
space charges are set up which interact with the rest of the electron-hole nairs 
along the n-track located within the originall) field fiee base m such a manner 


; This is so because charge collection by drift ceases when the slower carriers 
have been oicked up Only then will diffusion become predominant 


’kk ) — g 

; In good Si devices, recombination times are > 10 sec The influence of re- 
combination centers on charge collection may therefore be neglected as far as 
prompt charges are concerned. It cannot be neglected for the collection of 


delayed charges. 


t 


At the 


as to set up a strong electric field along the i-track within the base 
same time the original field of the depletion layer collapses Within a verv 
short time the electric field of the depletion layer lias moved into the base 
along the i-track, a ’'funnel" has opened up as it were, hence the name "field 
funnelling". Carriers within the base are therefore collected by drift rather 
than diffusion, a fact which was not anticipated before Earlier work on SEE 
emphasized diffusion as the main collection mechanism^^ Once the prompt 
charges have been collected (which takes about 10 ^ sec) the anginal field of the 
depletion laver will be restored and diffusion becomes gradual 1\ the more important 
collection mechanism field funnelling, being an essentially nonlinear space 
charge effect, comes about chen two conditions are ihjl (1) the original plasma 
column generated b\ the eneigciic particle must be of high densit\ (\ £ \ , and 
(2) part of the column must ti averse the depletion laser lhese at e the prei equ ls L te 
for the occurrence of field funnel.li.ng or enhanced prompt charge collection 
The funnelling may be qua 1 l tat ivc L\ described b' eomimng the one-dimensional 
junction case governed by equations (15) through (21) FLg la shows the 
energ>-band diagram tot a onc-d imens lonal n-n junction (t ~ 0-) An electric 
tield given by equation (L7) is present in the space-charge region Uc recall from 
semiconduc toi junction theory that the electriL fieLd in the fixed space-charge 
region establishes equilibrium between diffusive and drift tunent flow across the 
junction lor the situation m fig 1(a) where the applied bias potential is zero 
(V = 0) , the net current flow is also zero With no current flow there is no ohmic 
drop in the p and n regions ind therefore the field is zero outside the space-charge 
region. Tf* now, an ion-track traverses the junction, an el ect ron-hole plasma 
(^.p = v v n) is superimposed on the picture in lig la at t = 0 The electron-hole 
plasma is instantaneously affected by the electric field m the original space- 
charge region In this region, the plasma holes and electrons are accelerated 


to the left and to the right, respectively, by the space-charge field This 
initiates the ant l-SMnm^t r ic injection of holes and electrons into the p and n 
unions, respectively This sets up a double-injection plasma distribution, 
emanating from the space-charge region Ihe initial flow of carriers will serve 
to neutralize the space-charge, thereby cancelling the original field by some 
time t«0+ (Fig 1(b)) However, the excess plasma charge injected into the adja- 
cent regions will establish an electric field where it was zero before Tor 
Liu zero-bias i isc the fields established will be commensurate with Poisson^ 
cfju.it ion ^ ccj (8)) as deputed m \ ig 1(b) Even with zero-bias, a current 
wiLL 1 Low through lIiu contacts of the device in the opposite direction of positive- 
ly is cut rent l low (i e , is in i soLar tell) q he field distribution undergoes 
! either ciMiv es as time progresses, during which time the charge created bv the ion 
ir ilk is i ollci ted m the fonn of curient flow in the device The total collection 
oi ionized tlui f c is implemented b\ the combined effects of drift ind diffusion 
Houo'ti , the c tensior of the eltdric iield into the otherv iso field-ftee regions 
oi the device (it , what is termed "funnelling") bv the initial action of the 
original t leld in the spice-chirgc region on the ion-produced plasma, series to 
accCierite m t Lire tnc collection of the charge created b\ the ion track 

*lhe question mturilLx uises, "whv lias the effect of field funnelling no^ been 

discovered much earlier 1 " After iLL, p-n junction devices have been used routmeh 

since the c.irl\ soties as high energy particle detectors^. However the structure 

of these devices is such that the energetic particle track is confined entirely to 

the depletion Lner which has been made deliberately wide bv using high resistivity 

material (a low doping concentration) and a large reverse bias As can be seen 

13 -3 

from eq (15) a 1300 -cm sample corresponding to = 10 cm and a bias of 
300 volts gives a depletion layer width of W = 190 pm, rather large compared to 
the 2 m o r IC devices. At the same time the field free region of the base 


11 



is kept small, the back electrode being pl',ed right next to the depletion layer. 

Nonlinear space rhaige effects m particle detectors, somewhat akin to funnelling 

(here called "plasma effects") have nevertheless been observed and investigated 

17) 

b\ a number of authors . The plasma effect mav be described is follows The 

almost instantaneously generated electron-hole plasma along the track will be 

polari/cd at the outer edges of the plasma column owing to the electric field 

present m Lhe depletion layer Electrons and holes will be separated in a 

mrrow sheet and a nonzero space charge will be set up The space charge 

in turn creates an electric field inside the column which essentially screens 

out the 01 Lginal field lor the field configuration of eq (17), holes will 

escape at the end of the track and rove toward the cathode and a steady current 

will flow This process will continue until the whole plasma column is eroded 

Ihc plasm i time the tire it takes for the whole plasma to decay, can be 

computed with this model fSeibt at al and turns out to be 

1/3 


= L 32 ' 10 


-10 (x o V 


2VJ 

R 


[ sec ] 


( 22 ) 


m 

where 1 is given lv r cq (18) (in V/cm) , W by eq. (15) (in cm), N by eq (1) 
m U 

(m on S, K is the range of the plasma-generating f ast particle (m cm) and E^ 
is the lnitiaL energy of the incident particle (in eV) The model assumed 
a constant along die track. Ihis limitation lias been removed by Finch^\ 

It is here not the place to delve into all these considerations and analyses; 
suffice it to say that these attempts at dealing with nonlinear space charge 
effects suffer from serious drawbacks. The assumption of steady state erosion 
and the neglect of field distortion are but two examples of somewhat dubious 
simplifications made Co obtain an analytically tractable model. In fact, 
m searching the literature, we have found that there are scarce examples of 
analytically amenable problems dealing with the full nonlinearity of space 
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Fig 1(a) Band diagram for an equilibrium p-n junction 

in a one-dimensional configuration (V=0, t=0-) 



Fig. 1(b) Band diagram for npn-equil ibrium p-n junction after 
strike by single high energy ion (V=0, t=0+) . 



charge effects in a time-dependent setting All problems of this kind which 

have been solved analytically have exclusively dealt with single injection, 

18 ) 

i e , considering only the motion of one kind of carrier But field funnelling 

and related effects are basically double injection phenomena Electrons and holes 
must be treated on an equal footing No analytical methods are known to this 
author which can handle ncnlmear double injection phenomena therefore, computer 
programs have been devised to handle these situations One such piogram, developed 
by the IBM group , has already been mentioned^ Another code, developed b\ SRA^\ 
has also successfully tackled the field funnelling problem The reader is referred 
to the literature for details'*^ Like all computer codes dealing with complex 

nonlinear mathematical relationships, those devised recently for dealing with 
SEU’s also have their drawbacks They ate slow and expensive and it is therefore 
difficult to extract information in parametric form. One would like fer 
instance to know the time at whicn the maximum current induced by an *-stnke 
occurs, its magnitude, etc , as a function of doping concentration, applied volt- 
ages, etc. To extract this type of information from a computer code is of course 
rather awkward Consequently, attempts have been made m the recent past to 
incorporate field funnelling in heuristic models which are conceptually simple and 
at the same time accurate enough for engineering purposes These models will be 
discussed m the next section 


4 Comparison of heuristic SFU modelling with computer ca 1 cul at ions Ever 

since field funnelling has been discovered by computer simulation^^ people have 

been trying to construct simple theories incorporating the funnelling effect in 

20 ) 21 ) 22 ) 

one way or another The theories by Hu , McLean et al and Messenger 

are notable examples Let us firsL discuss the theory by Hu assumes that 

the funnelling effect produces a depletion region reaching far into the base 


/ 
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although the original depletion layer has collapsed Actually, the computer 
simulations of both Hsieh et al and Grubm et al^^ show that, although 

drift m the base constitutes the majoi means of charge collection for prompt 
charges, the base is far from depleted owing to the large residual carrier 
concentration within the energetic particle track during prompt charge collection 
Nevertheless, assuming a depleted region of width 


W = 


i + u 

n P 


Q(t)/qN 


0 * 


(23) 


where 


Q(t) 


( 


I(t) dr , 


(24 ) 


is the prompt charge coLLectcd at time t and is given by eq (I), Hu assumes 
that a potential difference develops, given by 


V(t) - \' b + V - q \ A V(t ) /2c r 0 , 


(25) 


f: 

. i 
\* 


where we recognize the last term of eq (25) to be the potential drop across a 
depletion layer valid for low level injection conditions. Eq (25) constitutes 
therefore a rather bold assumption Given the usual expression for the conductance 
C of a semiconductor, eqs (23) to (25) lead to an expression for the collected 
charge as a function of time 


Q(t) = Q h (1 - e" t/T H) 


(26) 


where 


and 


«H ■ ee 0 (v b + v > k a 1) 


, , 0 u + u 
1/2 p n 


N 0 ’ 


(27) 
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(28) 


*•) 

are tnc prompt charge and the collection time respectively Fig 2 shows a com- 
parison between eq (26) and the computer solution by Hsieh et al As can 

be seen, Hu’s simple model predicts a smaller collection time than the exact 
theor\ 


In their model for SEU’s McLean and Oldham do not address themselves to the time 

history of charge collection but proceed to derive an expression for the total 

21 ) 

collected prompt charge Ihey assume a funnel length L^ determined by the 

average drift velocity of electrons toward the junction and a charge ( ol lection 

time Ihe latter was estimated from lateral diffusion of the initial charge 

21 ) 

column and is given bv 


ML 


[9 


2 9 

o V 12 " D 



(V + Vl| 

b 


1/3 


(29) 


where L) is the imbipolar diffusion constant I lie total collected prompt charge 
was determined as (} = cj I and has the value 

/ 1/-^ c 1/6 i / i | i 1/9 

°M. * (3N n < V k + V > /8 ) (2 -n 9 /°) N a 7 (30) 

MI IJ b w n p A 

a rather unwieldy e>pression Using values for the various quantities occurring 

15 -3 

in eqs. (27) to (30) valid for our example (N =10 cm , 5 MeV u-particle strike, 
etc ) wt find 0^ =21 fC and = 65 fC for the prompt charges collected according 
to Hu’s model~^ and McLean-Oldham’ s model^^ , respectively. The charge collection 
times for the two models turn out to b^ = 180 psec and = 64 psec respec- 
tively. Incidentally, the total charge collected amounts to 100 fC for the 5 MeV 

*^In eq. (27) R signifies the range of the fast particle or the length of the 
device, whichever is smaller. 
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a-stnke and the collection time was 750 psec according to the computer calcula- 
tions of Hsieh et al As can be seen from the values for Q and T just quoted, 

for Hu's model, a third of the charge is collected m a time three times longer 
than for the model of McLean and Oldham But since both models are based on a li- 
ferent assumptions they do not necessarily agree. Furthermore, since the \ nder- 
lying assumptions, a fully developed depletion region m Hu's model on the one 
hand, and a funnel length determined by lateral diffusion in McLean - Oldham* s 

model on the other hand — assumptions which are not borne out by computer calcula- 
19) 

tions — are used, we must conclude that the two models discussed so far are of 

heuristic value at best. 


22 ) 

Yet another model has been proposed by Messenger He first proceeded to solve 

eq. (4a) for the case in which the fast ion penetrates the junction at normal 

incidence and is stopped within the depletion layer. He assumed that the column 

of electron-hole pairs is subjected to the full electric field of the junction 

Neglecting recombination, arguing that recombination times are much longer than 
*) 

transit times , he showed that the excess carrier density decays as 


6p 


-(]J E /w)t 
^ e ra 


(31) 


where p is a suitably averaged hole mobility, taking into account the dependence 


of p on E . E and W are given by eqs. (18) and (19) respectively. Another Lime 
p mm ' 

constant B ^ is then introduced ad hoc to describe the thermalization of hot elec- 
trons In the end, the transient current generated by the fast ion is shown to 
be proportional to the difference of two exponentials, viz , 


/ -(♦ E /w)t \ 

T I 171 “ 

I Me - e J 


( 32 ) 


*> 


This has been confirmed by exact computer calculations 


19) 
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Since Che mobility decreases with increasing field strength, we may use 
— 2 , 2 

i ~ 300 cm /V sec rather than the value of 600 cnG/V sec used earlier The 
characteristic time constant, eq (20), now becomes 


V = W/u E » 2 x 10 sec 
r m 


for our e%ample Defining t/ = \, the relationship (32) goes over into 


<(e*' - e"‘j= F( 


W’th = " 3he function F(\) exhibits the characteristics of the transient 

r 

current Qualitatively The current, initially being zero, rises fast towaid a 

maximum and than decreases, to become zero again, just as the computer simulations 

4 „ 15v , , . 19) . , , „ ,15) 

ot Hsien ) and ( n bin predict According to Hsieh , the current ma\imum 

15 ~3 

occurs for our example (8 V reverse bias, =10 cm , etc) at t= t =* 0 06 psec, 
corresponding to y - = 3. The maximum of F(x), eq (34), occurs at 

x = 05) 

m i-l 

Now, * must be greater than one m order for K(x) to be positive. Also, t > 1 

implies thermal izat ion times less than - , which is reasonable on physical grounds. 

However, with x = 3, eq (35) has no solution for < v 1 The root for this rncon- 
m — 

sistenc> lies m the neglect of nonlinear effects which are of prime importance 
here, as repeatedly stressed in this report The depletion la>er collapses m the 
neighborhood of the i-track, electric fields are generated m the formerly field - 
free bulk by space charges developed during the collapse of the origmallv strong 
built-in junction field This sequence of events, being caused bv nonlinear 
interactions between electrons and holes of the dense plasma created by ionization, 
cannot be modelled easily with the aid of simple pictures. After all, the 


funnelling effect has b en discovered by computer analyses, taking into account 
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the full nonlinearity of the system, which analytic methods were not then and are 
not now able to do 

5 Cone lusion An energetic particle traversing a p-n junction leaves in its 
wake a high density electron-hole plasma The subsequent decay of this plasma 
due to drift and diffusion leads to charge collection at pertinent electrodes 
which may upset the circuitry unduly. Since such events happen randomly and 
frequently given an adverse space environment, the necessity for understanding 
these events is obvious The basic equations governing the fate of the plasma 
(section 2) have been attacked for the problem at hand both by numerical com- 
putations and heuristic (analytical) methods Numerical methods are relatively 
slow and expensive The need for simple models which capture tne gist of the 
pertinent numerical work is felt quite vividly Owing to the nonlinearity of 
the formalism describing double injection phenomena (electrons and holes) as 
discussed in section 3, ab initio analytical methods cannot be devised success- 
ful] v Models, which have been discussed m section 4, have been shown to be 
inadequate This is due to excessive oversimplification, without which however, 
the models would become quite cumbersome, defeating their purpose The conclusion 
is inevitable Only numerical computations give reliable information on the 
time history of currents and collected charges generated by an energetic particle 
Only numerical computations take properly into account the nonhnearities which 
determine the shape of a current pulse. 
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